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Abstract. We define a quantity, the so-called purity fidelity, which measures the rate 
of dynamical irreversibility due to decoherence, observed e.g in echo experiments, in 
the presence of an arbitrary small perturbation of the total (system -I- environment) 
Hamiltonian. We derive a linear response formula for the purity fidelity in terms 
of integrated time correlation functions of the perturbation. Our relation predicts, 
similarly to the case of fidelity decay, faster decay of purity fidelity the slower decay 
of time correlations is. In particular, we find exponential decay in quantum mixing 
regime and faster, initially quadratic and later typically gaussian decay in the regime 
of non-ergodic, e.g. integrable quantum dynamics. We illustrate our approach by an 
analytical calculation and numerical experiments in the Ising spin 1/2 chain kicked 
with tilted homogeneous magnetic field where part of the chain is interpreted as a 
system under observation and part as an environment. 
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1. Introduction 

The relation between the rate of decoherence and the nature of dynamics is considered 
to be quite established JI], 0, |^. In particular, it has been shown that von Neuman 
entropy of the reduced (environment averaged) density matrix of an open quantum 
system (system -|- environment) grows (on a short time scale) with the rate given by 
the Lyapunov exponents of the corresponding classical dynamics. 

However this is expected to be true only under two rather severe conditions: (1) the 
system has to be initially prepared in a "coherent" state (i.e. minimum uncertainty wave 
packet) and (2) one has to be deep in the semi-classical regime of very small effective 
Planck constant h since the time-scale of expected quantum-classical correspondence 
of the entropy growth, namely the Ehrenfest time, scales as oc log{l/h). Indeed, 
performing real or even numerical experiments within this regime appears to be very 
difficult, see e.g. the numerical experiment in an iV— atom Jaynes-Cummings model @], 
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where a difference of the entropy growth between integrable and classically chaotic cases 
has been observed, but its quahtative nature does not become clear. 

On the other hand, by rejecting the assumption (1) on the coherent initial state 
but considering a random pure initial state instead, one can use orthogonal/unitary 
invariance and define simple random matrix models to analyze the time evolution of 
decoherence 0, and no qualitative and minimal quantitative differences between regular 
and chaotic dynamics were observed. Furthermore, we should stress that in quantum 
information science considering random initial states is potentially much more useful. 
Indeed, in order to use a massive parallelism of quantum computation one has to prepare 
the quantum computer system initially in a coherent superposition (pure state) of as 
many as possible elementary qubit (basis) states, but to contain maximal information 
this state can and should be any possible state and thus will behave more like a random 
state than like a gaussian wave packet. 

Based on a proposition by Peres 0, there has been much interest in viewing 
reversibility of a process with an imperfection on the reversed time evolution rather 
than on the state as discussed by Casati a/ 0. This allows to avoid some of the 
essential implications of linearity of quantum mechanics that trivialize the latter case. 
The related state correlation function is usually called fidelity. Its behaviour for chaotic 
systems on the Ehrenfest time scale has first been discussed in with similar findings as 
in ^, ^ . On the other hand it has been found recently [^, [11|] that on more relevant 
time scales related to the decay of correlation functions of quantum observables in the 
mixing case the fidelity decay is indeed exponential, but with a very different exponent, 
which is determined basically by the strength of the perturbation. This result was 
obtained for any state and, what is more, for the integrable case a faster gaussian decay 
was found in the thermodynamic limit ||^. More strictly speaking the result implied for 
small times or small perturbation strengths linear decay behaviour in the mixing case 
and quadratic decay in the integrable or more generally non-ergodic case, where 'small 
times' may be still large compared to the perturbative regime where all dependencies 
are quadratic anyway. 

There is a widespread sense that decoherence should follow correlation functions 



I4 , and we therefore expect that in the present situation decoherence should behave 



along the same lines as fidelity. Yet it is worthwhile to test this, as understanding 
decoherence is essential for the quantum information applications. The central issue of 
this paper is thus to test the evolution of decoherence of echoes (some aspects of which 
are relevant in the real spin-echo experiments |13|). Indeed it is crucial to determine if 
fidelity is a reliable measure. This is all the more true because of the counter-intuitive 
result that fidelity is higher for mixing systems than for integrable ones. Due to the 
great interest in quantum information problems we shall mainly focus on random states 
for the reason pointed out above. 

Basically we follow previous work on fidelity |TD|, ^ and express our results 



on decoherence by time correlation functions of quantum observables, though we shall 
have to generalize the concept of correlation function slightly. We shall use the purity 
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(which is one minus the hnear entropy or idempotency defect used in [Q) instead of the 
von Neumann entropy as a measure of decoherence, because purity, being an analytic 
function of the density matrix, is much easier to handle. Rather than introducing non- 
unitary time-evolution we follow the simpler way used in 0, ^ of partial tracing over the 
environment in a product Hilbert space after unitary time evolution. The kicked spin 
chain will serve as a specific example for our considerations, where analytic calculations 
for the integrable case can be carried further and numerics will confirm the validity of 
our approximations. 

The unitary propagator U of the total system (central system + environment) 
can be either a short time-propagator U = exp{—iHAt/h), or a Floquet map U = 
Texp(— z Jq dTH{T)/h) for a periodically time-dependent Hamiltonian H {H{t + p) = 
H{t)), or any abstract quantum map. We only assume that the total Hilbert space Ti 
(the domain of U) can be written as a direct product of two parts 

n = n,® He. (1) 

In what follows, the subscripts 'c' and 'e' will denote quantities referring to the central 
system and the environment, respectively, and to the total system when no subscript 
will be attached. 

As mentioned above we shall not introduce a non-unitary time evolution, but just 
consider the decoherence induced in the central system by its entanglement with the 
environment and then perform partial traces over the latter. In other words we shall 
test the stability of the disentanglement caused by inverse time-evolution operator if this 
operator is perturbed. This perturbation can generally be described by some self-adjoint 
operator A over Ti and some perturbation strength S to yield 

Us = Uexp{-iA6/h). (2) 

Our perturbation is static (the generator A has no explicit time-dependence) and cannot 
be associated to noise since it is not due to the coupling to the environment (which 
is explicitly included from the outset) but stems from the imperfect (or unknown) 
description of the total hamiltonian. 

First, let us consider coherent quantum evolution of the total system, and define 
the fidelity as an overlap between states evolving under unperturbed and perturbed 
time-evolution 

Fit) = {Mt)\m) = mtu'i^p) (3) 

where our time t is an integer. Fidelity has become a standard measure of instability of 
quantum computation |T^, but has also been considered in a more abstract context as 
a measure for the instability of quantum dynamics [|, ||, ^ |10|, |T6|, 0, |lT], |19|, p^, [20|]. 
It has recently been pointed out ^ that fidehty is intimately related to the decay 

of correlations, and that, surprisingly enough, fidelity decay is faster the slower decay of 
correlations is. However, fidelity is a property of a pure state of the total system which 
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is typically not accessible. What can defacto be measured is only information relating 
to the reduced density operator 

p,it) = tr,\m){m\ (4) 

where tre is a partial trace over the environmental degrees of freedom Tie- 

We now have to obtain an extension of the concept of fidelity in order to measure 
coherence properties of the reduced density matrix. Assume we prepare our system in 
a pure product (disentangled) state 

|^) = |^c)®|V^e), (5) 

such that the reduced density matrix is also pure Pc(0) = \4'c){4'c\- Then we propagate 
our system for some time t, after which the state U'^lip) becomes generally an entangled 
superposition of system and environment states. We then invert time, i.e. we change 
the sign of the Hamiltonian with a small inaccuracy described by the operator A (H), 
and propagate the system backwards for the same amount of time arriving at the final 
'echo' state 

\m = w'u'i^). (6) 

The overlap between the initial and the final state would be just a fidelity 
F{t) = {(j){t)\'ip) , however here we are interested only to what extent has the final state 
disentangled from the environment i.e. how much the final reduced density operator 

pr'°W = tre|0(t))(0(t)|. (7) 

deviates from a pure state. This is best quantified in terms of a purity, trp^, which 
is equal to 1 for a pure state and less than 1 otherwise. The minimal value of purity 
is l/N'c where is the dimension of Tic] note that this limiting value can only be 
reached if the dimension of the environment AC tends to infinity. For finite dimensional 
environment space we find later a more accurate limiting value ( |33D (see also f^). We 
shall therefore study, the purity fidelity defined as the purity of an echoed reduced state 



Fp(t) = tr, [pT'^^it)]' = tr, [tre {Ui'U'\^,)\tP,){tP,\{^,\U-'U's)'_ 

In section 2 we shall make our main theoretical predictions on purity fidelity in 
relation to ergodic properties of dynamics and in particular to the correlation decay. 
In section 3 we shall apply our results in the quantum spin chain model, namely the 
Ising spin 1/2 chain kicked with a tilted homogeneous magnetic field, which exhibits all 
qualitatively different regimes of quantum dynamics ranging from integrable to mixing. 
In section 4 we shall discuss some implications of our results and conclude. 

2. The relation between purity fidelity and correlation decay 
2.1. Linear response 



2 



Following |Ty, |TT] we start from time-dependent perturbation theory (linear response) 
and expand the purity fidelity in a power series in the perturbation strength. It is 
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convenient to write the complete basis of 7i as \ j, u), where Latin/Greek indices running 
over J^cc = dim7ic,c values denote system/environmental degrees of freedom, such that 
the initial state is always designated as |1,1) = li'c) ® li^e)- First, let us rewrite the 
purity fidelity 

Fp{t) = tr, [tre [Mt\^,)\^e){^c\{ipe\M})f = (9) 
= J2 (j,/i|Mi|l,l)(l,l|M/|fc,/x)(fc,z/|M,|l,l)(l,l|M/|j,z/) (10) 

in terms of a unitary fidelity operator 

Mt = Ui'U'. (11) 

Second, we observe @ that the fidelity operator can be rewritten in terms of the 
perturbing operator in the Heisenberg picture At = U~'^AU^, namely 

Mt = exp{iAoS/h) exp{iAi6/h) ■ ■ ■ exp{iAt-i6 /h). (12) 

This expression can be formally expanded into a power series 

oo -rnxm t 

^* = l+Ezni^^E At,At,---At^, (13) 

m=l ti,...,t™=l 

which always converges (in the strong limit sense) provided the generator A is a bounded 
operator. The symbol T denotes a left-to-right time ordering of operator products. 
Third, we truncate the expression (|T3D at the second order 5^ and plug it into expression 
(y). After a tedious but straightforward calculation we find 

ox2 t-1 jj^lu^l 

FAt) = i-^ E EE(i,i|A'|j,^)(j,HA"|i,i) + (^('5')- (14) 

t',t"=0 j V 

The RHS of this linear response formula indeed looks much like a time correlation 
function, however with funny exclusion rules on the state summation. It can be written 
more elegantly in a basis independent way as 

Fp(t) = 1 - ^tr [pSAtpSA] + 0(5^), (15) 

where = EtrJoA', P = \^){^\, and p = (1^ - |Vc)(^c|) ® (le - |V'e)(^e|). 

Our expectation is that the properties of purity fidelity decay will mainly depend on 
dynamics, i.e. the behavior of the operator SAj, and less on the detailed structure of the 
initial state encoded in p and p. Let us first discuss the limiting qualitatively different 
cases of dynamics: 

2.1.1. Regime of ergodicity and mixing. In the regime of ergodic and mixing quantum 
dynamics |^T|, the reduced transport coefficient 

^= liniltr[pSApSA] (16) 

can be estimated in terms of the usual Kubo transport coefficient 

^ t-1 

^=1™- tr[pA'A"] (17) 

t't"=0 
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which is obviously finite cr < oo in the case of mixing dynamics and sufficiently strong 
decay of time correlations. Namely one can easily prove that < a < a since all 



the terms of a missing in a, as well as all the terms of a when expanding along (H), 
are non-negative. For times t larger than a certain mixing time t^i^, t > tmix, i-e. a 
characteristic time scale on which the limiting process (|T6|) converges, we thus find a 
linear decay of purity fidelity 

Fp(t) = 1 - ^at + 0(6') = 1 - + 0(6'), (18) 

iT' Tern 

on a time scale 

rem = ^oc 5~\ (19) 

2.1.2. Non-ergodic regime. In the opposite regime of non-ergodic (e.g. integrable) 
quantum dynamics, the non-trivial (e.g. different than a multiple of identity) time- 
averaged operator exists 

A = lim -^At, (20) 

so that for times larger than a certain averaging time t > tave, i-e. a characteristic time 
scale on which the limiting process ( pOD converges, we find a quadratic decay of purity 
fidelity 

Fp(t) = 1 - -^Dt^ + 0(6^) = l-2(^) + 0(5^), (21) 

with a coefficient which we name as reduced stiffness 

D = ti[pApA\ (22) 
on a time scale 

f,e = oc 6-\ (23) 

We note that our reduced stiffness coefficient is again rigorously bounded by the common 
(non-reduced) stiffness D = tr[pA^], namely D < D. 



2.2. Beyond the linear response 

Next we extend our linear-response results to the regime, where the value of purity 
fidelity becomes considerably lower than 1. Again we consider the two qualitatively 
different cases of dynamics: 



2.2.1. Regime of ergodicity and mixing. Here we assume, in full analogy with the 
derivation of fidelity decay [TT|] , that also the property of quantum n-mixing holds (in 
the asymptotics A/" oo), namely that all the higher n— point time correlation function 
factorize so that the 5— expansion of the fidelity operator can be summed up and in the 
weak limit sense we obtain 

e^/^^-^Mt ^1, as t > tmix, where Tem = (24) 
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This statement is equivalent to the statement shown in |TT] that for arbitrary (pure 
or mixed) state p, fidelity decay F{t ^ tmix) = ti pMt = exp(— t/rem) is independent of 
the state. In case of a semi-classical situation of small effective value of h the limit (^) 
starts to build up only after the Ehrenfest time ^e, = ln(l/^)/A, where A is the 
largest classical Lyapunov exponent, so then t^i^ = ^e- Of course, (0) never holds in 
the strong limit sense as the fidelity operator Mf is unitary. Plugging ( P^ into formula 
(I) we obtain for the asymptotic exponential decay of purity fidelity 

Fp(t) = exp(-4t/rcm), as t > tmix- (25) 

Comparing eqs. (^) and ([T8|) in the asymptotic regime t — >■ oo while approaching 
5 —> so as to remain in the regime of linear response we find that we must have 

strictly fem = 1"em, 1.6 

a = a. (26) 

Of course, this argument is correct only for sufficiently strong mixing, e. g. for 
exponential decay of correlations. 



2.2.2. Non-ergodic regime. Again, in analogy to the simple fidelity decay in non-ergodic 
situations^, we can rewrite the fidelity operator in terms of the time averaged 
operator A, 

Mt exp{iAt6/h), as t > tave- (27) 

Plugging the fidelity operator (^) into formula (^ we find that, as a manifestation of 
non-ergodicity, the behaviour of purity fidelity generally depends on the structure of the 
initial state. 



2. 3. State averaged purity fidelity 

As we stated in the introduction, we are interested in the behaviour of a random initial 
state of the central system and as far as the environment is concerned we as well average 
over states to refiect our ignorance of the latter. Let us therefore start with a random 
product initial state (|) implying that the states \ipc) and \ipe) are random A/'c and Afe 
dimensional vectors, whose components (jlipc) and (yulV^e) are, in the limits A/'c oo 
and A/'c —>■ oo, independent complex random gaussian variables with variance l/A/'c and 
l/A/'e, respectively. Denoting averaging over random product initial state as (.)^ one 
can easily average the linear response formula ([l^ 

expressing the purity fidelity decay in terms of a sum of reduced correlation function 

Cp(t', t") = {AeA„) + - ((A')e(A")c)c - ((A')c(A")c)e, (29) 

where (.)c = (l/K) tic (.), (.)c = (l/K) tr, (.), and (.) = ((.)c)c = (l/AT) tr(.) 
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Again, in ergodic and mixing regime, where correlation functions (^) decay fast as 
\t" — 1'\ ^ oo, the average purity fidelity exhibits initial linear decay ( p!8D with average 
transport coefficient 

(a)^= hmi Cp{t',t"). (30) 

~* t',t"=0 

Since, for sufficiently strong mixing we have (p6[), (cr)^ = limt^oo ^ Ylt~t"=o^i'^' ~ '^") — 
{a)^, where C{t) = {A^A) - {A)', then lim,^„, i El-)'=o{C{t' - t") - Cp{t', t")) = 0. 
We note again that all the results referring to the ergodic and mixing regime implicitly 
assume the limits J\fc,Me — >■ oo to be considered before t — oo. 

On the other hand, in non-ergodic regime we have quadratic initial decay (|2T| ) with 
state-averaged reduced stiffness 

{Dh = +T)(!v; + 1) (^^'^ + " ^^^^'^^ " ^^^^'^^) • ^^^^ 

Without making reference to any of the two extreme cases the general expression for 
the purity fidelity (^) may be simply state averaged while keeping the fidelity operator 
completely general: 

^ "^ '^'^ ~ {K + im + 1) ^ KiK + miK + 1) ^ ' 

Y ( {P, a\Mt\j, iy){j, /i|M/|g, (3\Mt\k, z/|M/b, (3) 

j,k,p,q,a,l3,tj.,u 

+ {p, a\Mt\j, u){k, u\M}\q, a){q, P\Mt\k, /i|M/|p, P)) 

From ([3^ ) one observes that (averaged) purity fidelity does not decay to zero. It 
has a finite saturation plateau value F*, which can be estimated by assuming that 
after a long time t and for sufficiently strong perturbation 6 (i.e. such that the 
eigenstates of Us look random in the basis of eigenstates of U) the fidelity operator 
Mt becomes a matrix of independent gaussian random variables. Then the average of 
the product of four matrix elements ( |52D can be estimated by a pair-contraction rule 
with ((p,a|Mi|j,z/)(g,/3|M/|fc,/i)) 

p. ^ K + K + ^ + Oji/K) + Oji/K) . . 

(A4 + i)(A4 + i) • ^^^^ 

We have defined a reasonable measure for decoherence for the echo process and 
have determined that purity decays faster for near integrable system than for a mixing 
one. This fact is quite surprising in view of previous results on decoherence, but not 
so much if we consider that it follows closely the results obtained for fidelity. We will 
proceed to apply our findings to a kicked Ising spin chain. 



3. Kicked Ising chain 



We shall now consider application of purity fidelity decay in a class of model systems 
namely one-dimensional spin 1/2 chains. We consider a ring of L spins described by 
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Pauli operators cr", j G {0, 1, ... L — 1}, a G {x, y, z}, with periodic boundary conditions 
0""^^ = (y'j- In particular we concentrate on the example of Kicked Ising (KI) model 
with the hamiltonian 

^Ki(t) = E [j^^^^Ui + ^pit){h.a'; + /i,(T])} (34) 
i=o 

where 5p(t) = Ylm=-oo S(t — mp) is a periodic delta function, generating the Floquet-map 
U = exp ^-iJ, E S'^i+i j exp l^-i E(^xfTj + /izfTp j , (35) 

where we take units such that p = h = 1, depending on a triple of independent 
parameters {Jz.hy^.h^)- KI is completely integrable for longitudinal {h^ = 0) and 



transverse {h^ = 0) fields ||22[, and has finite parameter regions of ergodic and non- 
ergodic behaviors in the thermodynamic limit L ^ oo for a tilted field (see fig. p. The 
non-trivial integrability of a transverse kicking field, which somehow inherits the solvable 



dynamics of its well-known autonomous version |23], is quite remarkable since it was 



shown that the Heisenberg dynamics can be calculated explicitly for observables 
which are bilinear in Fermi operators Cj = (aj — iaj) 11/^"' with time correlations 
C{t) decaying to the non-ergodic stationary values D as |C(t) — D| ~ In order 

to test our predictions of purity fidelity decay by explicit calculation and/or numerical 
experiment we consider a line in 3d parameter space, with fixed = l,h^ = 1.4 and 
varying exhibiting all different types of dynamics: (a) integrable for = 0, (b) 
intermediate (non-integrable and non-ergodic) for = 0.4, and (c) ergodic and mixing 
for hz = 1.4. In all cases we fix the operator A = M := J2j <^J which generates the 
following parametric perturbation of KI model with 

.^^.^^iJAlJpS^S^OiS'), (36) 

^ ''--''''-^ ^J"'°^'''' i + 0{y-), where h = ^hfThi 

and vary the size L and the perturbation strength 6. 

For the sake of mathematical simplicity of the model we do not want to introduce 
extra degrees of freedom to simulate the environment. Instead we logically split 
the spin-chain by assigning a subset JT^ of spins to the central system and the 
complement = Jc oi = L — spins to the environment. Then we have 

our numerical experiments and explicit analytical 
calculations we shall consider three special cases (in the first two cases we have to 
assume that L is even): 

• (A) Alternating subchains, where every second spin is assigned to the central system 
J^^ = {0, 2, 4 . . . , L — 2}. This situation is supposed to be a good model of a 
situation where all central system's degrees of freedom are directly coupled to the 
environment. 
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• (C) Connected subchains, where a connected half of the chain is assigned to the 
central system J'^ = {0, l,2...,L/2 — 1}. Here the central system is coupled to 
the environment just at two ending points, namely j = and j = L/2 — 1. 

• (S) Single spin, where the central system consists of a single spin J'^ = {0}, i.e. a 
two- level quantum system coupled to a correlated many-body environment. 

Our explicit calculations will be performed for a random initial state, in other words 
only state averaged purity fidelity will be computed. 



3.1. Exact calculation for the integrable case 



Here we show how to extend our analytical approach to solve purity fidelity decay 



in the integrable case h^, = 0. We start by a finite size version of the formalism, however 
the final results will be most elegant and simple in the thermodynamic limit L ^ oo. 
The main quantity which we want to calculate is a state-averaged reduced stiffness {D)^ 
(Pip. As a first step in our calculation we have to construct a time average operator 
A associated to the perturbing operator A. This can be done, as described in [^, by 
means of the invariant space of the adjoint map U^'^, which is a unitary operator defined 
on a linear space of (bounded) observables as U^'^A = U'^AU. The complete basis of 
this invariant space can be constructed in terms of two sequences of 2L operators, Un, 
Vn,neZ2L, namely 



Un = -al n = 0, (37) 



where {o'J)k '■= IlfLi for k > 1, (cr|)o := 1, satisfying a Lie Algebra 

=0, (38) 

\Umi '^n\ '^'i'iUm+n ^m— n)- 

We note the orthogonality 

{UnU^) = (KKn) = {UnVm) = 0- (39) 

Now, the Floquet-map ( pSD can be written in terms of the elements of the algebra (p^) 
only, namely U = exp{—iJUi) exp(z/ixf/o); and the complete basis of an invariant space 











n 


>1, 










n 


= 0, 


j=0 


[ 






n 


<-l, 








n-lO"j-|-ru 


n 


>1, 


e| 








n 


= 0, 


j=0 






-n-lCrJ_„, 


n 


<-l, 



can be, for arbitrary finite L, written as "^(^i) with (fi = ^,lE 'Z2L, and 

2L-1 

^{if^)= $:(ix(^,)e^'^^'+^(-V'Oe-^"^0-^^n, (41) 



n=0 
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K [/_„,2-i/2(v;-v^„„)), 

u{(p) = (cot a - cot /3e~"^, cot a - cot /3e*'^, v^z sin^j), 

where we introduced the angles a = 2J, /5 = 2hx. Note that '^{(p2L-i) = ^(v^O and the 
basis (of non-equivalent vectors / = 0, 1, . . . , L) is orthogonal 

{^^M^i^i)) = 2L^\u{^k)\\6k,i + 6k,2L^i). (42) 

The time average of an observable can be simply constructed by means of a projection 

namely for the magnetization A = M = Uq v/e have 

In the limit L —>■ oo the phase variable becomes continuous ip G [0, 27r) and the sums 
(|^,|^) go over to integrals over ip. The first term of (|3l|) is the stiffness D = (M^) which 
can be calculated explicitly in the limit L ^ oo Finite size corrections are proven 
to be smaller than any power in 1/L, i.e. they have to be exponentially small, since the 
argument of the sum which is approximated by an integral over ip is an L-independent 
analytic function of e*"^ 

max{| cosal, | cos/3|} — cos^ P 
sin^ P 

In order to calculate the remaining terms of the reduced stiffness ( |31| ) we have to compute 
the reduced inner product of a pair of operators, namely 

{{A)j{B)j)j, = 2~^^'^-'^^hTj,{tTjAtTjB) (46) 

where J' is some subset of Z^, (either J'c or J'^) and J'' = \ JT". We start with the 
basis operators f/„, Vn, for which we find by direct inspection of the structure (|37D , that 

{{U^)j{Un)j)j' = {{V„,)j{Vn)j)j' = LfJ^^, (47) 
mn.)j{Vn)j)j' =Q. 

Here the structure function is defined by the number of sequential sets C(/i,/2) = 
{/i, /i + 1, . . . , I2} in Tjl of length \n\ + 1 which do not intersect the set J , divided by L 

fn = ^\{jeZL; C{j,j +n)nj = ^}\, (48) 
and is extended to Z2L by putting /__„ = /2L-n := fn- Using a form factor of sorts 

2L-1 

F(^) = Y: /ne*"^ (49) 

n=0 

we can directly write the reduced inner products within the invariant family (p|) 
generalizing (^2]), namely 

= 2Lu{^r ■ u{^) (Ficp - + + ^')) . (50) 



D = L ^ \l ^ / '^ + O (exp( -const L)) . (45) 
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The calculation of (pi]), the reduced inner product of a time-averaged observable A with 
itself, {{A)'j) jt is of central importance. It can be computed straightforwardly using the 
expansion (^31) and reduced inner products (|50|). For the magnetization M we obtain 



X 



x{F{^i-^k) + F{^i + ^k)). (51) 

The factor L has been taken out in the definition of the coefficient Gi;' in order to 
make it conveniently size-independent for the thermodynamic limit. In order to proceed 
with explicit calculations we have to derive the form factors for the three different cases 
(A,C,S). For each of them we have to compute two generally different reduced inner 
products with the form factors, Fc for J = J'c, J' = J'e, and Fe for J' = J7e, J7' = J7c, 
corresponding respectively to the last two terms of (pTF). 



• (A): Here, due to symmetry both structure functions are identical = |5n,o giving 
FH^) = F^^) = \. 

• (C): Again, we have a symmetry between the central system and the environment, 
so = max{i - ^,0} and F^{ip) = (v?) = -ndi/ii^) where 5m{^) is a fat 
periodic delta function 

U^) = ^{"^^^^]\ (52) 
27rm y sin(y9/2) j 

with the limit limm^oo ^m(v^) = Y^'kL-oo ~ ^nk). 

• (S): Here, we need to consider two different structure functions f^^ = (L — 1 — 
\n\)/L, = l5n,o giving for the form factors F^{(f) = F^{(f) = ^^^^|^5l_i((^). 

Due to linear dependence of G^? on F{ip) ( pT| ) we need to compute only two types of 
the reduced inner products besides the non- reduced one (H), namely for a constant 
form factor F{(p) = 1 and for a fat-delta form factor F{ip) = 6mi^)- Tedious but 
straightforward calculation gives (within accuracy 0(exp(— const L)) beyond all orders) 

2 cot P sin^(a/2); cos a > cos P, cos a + cos P > 0, 

— tan(/3/2); cos a < cos /3, cos a + cos /3 > 0, 

cot(/3/2); cos a > cos/5, cosa + cos/5 < 0, 

2 cot P cos^(a/2); cos a < cos P, cos a + cos P < 0, 

^'- = ^'-^' ^^^^ 
The state averaged purity fidelity in the regime of linear response (|2^) can be therefore 
written as 

(^^m-^-J^^^§^f^'^Oin (54) 



Gi 



G 
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where the coefficient in the three different cases of interest reads 
= 2710 s — Gi, 

= ^, (55) 
= ^ (2nGs - Ci + i 

It is very interesting to note that only in case (A) where all environmental degrees 
of freedom are coupled to the central system and vice versa, the resulting expression 
for purity decay coefficient G is non- vanishing in the thermodynamic limit L oo, 
so the purity-fidehty decays qualitatively in the same way as the (non-reduced) fidelity 
1^, namely the exponent of quadratic decay in time has the same L dependence. In 
the other two cases (C,S) where the interaction between the central system and the 
environment is local one has a subtle cancellation of leading order terms giving the 
resulting coefficient which vanishes as oc 1/L in the thermodynamic limit thus making 
the purity-fidelity to decay still quadratically but on much longer time scale (by a factor 
of L) as compared to the time scale of fidelity decay Pj. Of course, our analytical results 
for the integrable case only give us a time-scale of the purity-fidelity decay 



1 



(K+im+i) .... 

but they cannot tell us anything about the global behaviour of {Fp(t))^ beyond the 
regime of linear response. These analytical results are clearly confirmed by direct 
numerical simulations reported in next subsection. 



3.2. Numerical calculations in the general case 

In the generally non-integrable cases of kicked Ising model we calculate the partial 
correlation sums 

Sj{t) = (57) 
and the total purity correlator 

Mt) = j E Cp{t',t") = S^{t) + S^,{t)-Sj{t)-Sj>{t) (58) 

^ t',t"=0 

by means of numerical simulation. With these we can calculate purity fidelity and 
compare to the directly numerically simulated purity fidehty decay (and to the analytical 
calculation in the integrable case). Note that S(D{t) = -^([SA^]^) = j^Yft^t"=o{^^t'-t") 
is the usual (non-reduced) integrated correlation function and S'z^it) = due to 
(At) = {A) = 0. We work with a random initial states, hence in the linear response 
regime (sufficiently small 5) we have an identity (|2^) connecting the state-averaged 
purity fidehty decay to the purity correlator. Asymptotic (5 0) exactness of this 
relation has been carefully checked in all different cases of interest since it provided a 
crucial test of our numerical procedures. 
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In the first set of numerical calculations we check the decay of correlations, either 
with respect to the full trace (non-reduced) inner product or with respect to the 
reduced inner product, while in the second set of calculations we explicitly compute 
the purity fidehty decay {Fp(t))^ and compare the exponents with our predictions 
based on correlation decay. Calculations are performed with three different system 
sizes namely L = 12, 16, 20. Since we want the perturbation strength to have certain 
size L- independent effects we scale it by fixing 6' = 6^Jl/Lq when varying L where we 
choose Lq := 24. For example, in such a case the predicted exponents of the fidelity 
decay do not depend on L. 

In fig. ID we examine the decay of non-reduced correlations of the magnetization 
C(t) = {MMt)/L for the three cases: (a) integrable = 0.0, (b) intermediate 
hz = 0.4, and (c) ergodic and mixing hz = 1.4. We find that C{t) has a non- vanishing 
plateau (stiffness C = D/L) in the integrable and intermediate non-integrable case, 
where the correlation function in integrable case agrees excellently with the analytical 
result for L — > cxd while in ergodic and mixing case we find exponential decay of 



correlations. We note an interesting distinction between correlation decays in integrable 
and intermediate cases, namely in the integrable case (a) the relaxation of C{t) towards 
the plateau value D/L is a power law ~ t^^/^ whereas in the intermediate case (b) it 
looks like an exponential (see inset of fig. |I]b). 

In fig. 1^ we show partial correlation sums Sj{t) for the ergodic and mixing case 
(c) for different reducing sets JT" appearing in the three cases (A,C,S) of divisions. We 
note that Sj{t) increases linearly in t only for the non- reduced case J7 = where (non- 
reduced) correlation function is homogeneous in time, whereas in all other cases Sj{t) 
increases slower than linear, in fact in most cases it quickly saturates to a maximum 
value which does not increase with increasing size L. Therefore, for large system sizes 
L in the ergodic and mixing situation the purity correlator becomes determined by the 
total (non- reduced) correlator only [confirming eq. (p6|)], namely Sp(t) = st, where 
s = 2.54. This is illustrated in fig. ^. 

In figs. ^,|^ we show analogous results on partial correlation sums and purity 
correlators in the integrable case (a). Here, Sp{t) oc t"^. Numerical results are compared 
with analytical expressions (0,^) and the agreement is very good. In fig. |^ we also show 
the purity correlators in the intermediate regime (b) where the results are qualitatively 
very similar to the integrable case (a), namely we have the quadratic growth of purity 
correlators Sp{t) oc t^. 

In the second set of numerical experiments we calculate the purity fidelity Fp{t) 
with respect to a random initial state and average the result over a sufficient number 
of realizations of random initial state such that statistical fluctuations are negligible. We 
note that Fp{t) is self-averaging, namely by increasing the dimension of the Hilbert space 
i.e. increasing the size L the purity fidelity of a single random initial state converges to 
the state average {Fp{t))^. Here we are interested not only in the linear response regime 
but also in the global behaviour of purity fidelity, so we chose to display the results for 
the (scaled) perturbation strength 6' = 0.01. 
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Figure 1. Correlation decay for three cases of Kl: (a) integrable hz = 0, (b) 
intermediate = 0.4, and (c) ergodic = 1-4, for different sizes L = 20,16,12 
(solid-dotted connected curves, almost indistinguishable in (a,b)). Circles (a) show 
exact L = oo result. Chain lines are theoretical/suggested asymptotics (see text). 

According to our relation (p^) and the behaviour of purity correlators we predict 
(and find!) that purity fidelity decay is faster in integrable and intermediate cases 
(a,b) than in the ergodic and mixing case (c). This result is most clear-cut in the case of 
division (A). See fig. |^ for a comparison of purity fidelity decay in integrable and mixing 
cases. We have found very good agreement with the predicted exponential decay (p5|) 
in the mixing case (c), whereas in the integrable regime (a) numerical results suggest 
a global gaussian decay of purity fidelity (similarly to the fidelity decay 0) with an 
analytically computed exponent (0). Due to finite Hilbert space dimension we find 
a saturation of purity fidelity for very long times at the plateau value which has been 
computed in eq. (^5]). In order to avoid the trivial effect of a constant term F* in ( |5^ we 
subtract it and in the following plots show the quantity \Fp(t) — F*\/ {1 — F*) . In fig. ^ 
we show that purity fidelity decay in ergodic and mixing regime (c) is independent of the 
type of division (A,C,S). In figs. we show an analogous comparison of purity fidelity 
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Figure 2. Partial correlation sums in the mixing regime /i^ = 1.4 for different divisions 
TLl = J ^ J' ■ The structure (alternating A vs. continuous C) and/or the number 
of elements of the set J' is indicated in the label near each triple of curves (full for 
L = 20, dashed for L = 16 and dotted for L = 12). 




Figure 3. Purity correlator 5p(f) in the mixing regime = 1.4 for different cases 
(A,C,S, and total (non-reduced) correlator, for full, long-dashed, short-dashed, and 
dotted curves, respectively) as compared to the asymptotic linear increase (chain line). 
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Figure 4. Partial correlation sums in the integrable regime hz = 0.0 compared to the 
analytically computed coefficients for different divisions indicated by labels describing 
the structure (con. /alt.) and number of elements of the set J' . 




Figure 5. Purity correlator Sp{t) in the integrable regime — 0.0 compared with 
an analytical expressions (sampled symbols) for different cases of divisions (A,C,S). 
fn the inset we emphasize the quadratic increase by plotting S'p [t) j and comparing 
to theoretical coefficients (55) (dotted). 
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Figure 6. Same as in figure (^) but for tlie intermediate case (b). Since we 
have no analytical predictions here, we compare the data by best fitting quadratic 
functions (sampled symbols, and dotted in the inset) whose coefhcients are later used 
for comparison to purity fidelity decay. 

decay in integrable and intermediate cases (a,b) for different divisions and show that, 
in all cases, short time behaviour is well reproduced by the linear response coefficients 
given by purity correlators. For longer times a global gaussian behaviour with theoretical 
exponents work quite well (in particular for case (A)). 

4. Discussion 

We have analyzed the properties of decoherence in the framework of the question of 
reversibility with a perturbed time reversed Hamiltonian. Up to now only the state 
correlation function, usually called fidelity, has been discussed in this context. The 
calculations were carried out in the framework of partial tracing over the environment 
and unitary time-evolution. The remaining density operator on the subspace was 
analyzed in terms of the trace of its square, usually called purity. It is more convenient to 
use this quantity rather than entropy itself, because its analytic form allows to perform 
explicit calculations to a much larger extent than would be the case for the entropy, 
where the logarithm complicates things. 

We apply the techniques developed for the calculation of fidelity to a Hilbert space 
which is a product space. Keeping track of the respective indices we can readily perform 
the partial traces needed, and find that they relate to slightly modified correlation 
functions of quantum observables. This allows to compare the results for purity-fidelity 
with those of fidelity 0. 

As for the fidelity it becomes apparent that the decay time of correlations is a 
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Figure 7. Purity fidelity decay: comparision between mixing — 1.4 and integrable 
regime = 0, at scaled perturbation S' ~ 0.01. Sampled symbols give the predicted 
exponential (^5|) in the mixing regime and a suggested gaussian with a theoretically 
computed time-scale (56|) in the integrable regime, while horizontal chained lines give 
the plateau values (p3). 




Figure 8. Purity fidelity decay in mixing regime — 1.4 for different types of 
division and different sizes L (indicated in a legend) and 6' = 0.01. Theoretical decay 
(for L oo) is given by sampled symbols. 



Decoherence of spin echoes 



20 




50 100 150 200 



t 

Figure 9. Same as in fig. || but for integrable regime (a). Note that tlieoretical 
(gaussian extrapolated) curves for cases (C,S) are practically indistinguishable. 




relevant short time scale, and this induces for mixing systems an exponential decay 
with a time scale that is mainly related to the strength of the perturbation. This 
implies a linear decay after the end of the trivial short time 'perturbative' regime 
which leads to a quadratic decay always. For integrable systems on the other hand 
correlations do not decay and thus the quadratic decay of purity- fidelity survives beyond 
the perturbation regime. This amounts to the central result that decoherence will be 
faster for an integrable system, than for a chaotic one in strict analogy to the findings 
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for fidelity. 

The techniques developed are applied to kicked spin chains where analytic results 
for the integrable case and numerical results for the correlation functions of observables 
can easily be obtained. In the thermodynamic limit we obtain exponential decay of 
purity-fidelity for the mixing case and gaussian decay for the for the integrable case 
with uniform coupling to the environment, which is simulated by assigning alternatively 
one spin to the environment and one to the central system. If the spin chain is simply 
cut in two or if a single spin is associated to the central system we find deviations 
from the gaussian shape, but the decay rate associated with the quadratic behaviour 
at short but non-perturbational times is correctly reproduced thus maintaining the fact 
that even in these situations coherence decays faster in the integrable case than in the 
mixing one. 

The central message is thus that decoherence in echo situations follows fidelity. 
Indeed the reversibility of decoherence under perturbation of the reversed Hamiltonian 
is better for mixing systems than for integrable ones, and the time scales on which this 
occurs depend sensitively on the strength of this perturbation. 
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